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Abstract 

We model the transition of a system of quarks and gluons at high energies to a 
system of quarks and mesons at low energies in a consistent renormalization group 
approach. Flow equations interpolate between the physics of the high-temperature 
degrees of freedom and the low-temperature dynamics at a scale of 1 GeV. We also 
discuss the dependence of the equation of state on baryon density and compare our 
results with recent lattice gauge simulations. 



The equation of state of QCD in extreme conditions of high temperature 
and/or high baryon density is at the focus of interest in heavy- ion collisions 
and astrophysics, i.e. neutron stars and supernovae. The current heavy-ion ex- 
periments at SPS (CERN) and RHIC (BNL) have challenged theoretical work 
on strong interactions at finite temperatures, due to the possible observation 
of new phases of matter. In particular, one expects a phase transition which 
separates a low-temperature region of hadrons from a high-temperature region 
of quarks and gluons, where both confinement is absent and chiral symmetry 
is restored. During the last years various direct numerical lattice simulations 
of QCD at finite temperature and finite baryon density have been performed 
[1-3]. These computations are still limited to rather small lattices, therefore 
an alternative approach using an effective field theory is useful. Previously, 
we have discussed a linear a-model with two massless quark flavours [5,6] 
which is treated with the renormalization group to capture the dynamics of 
the long-range fluctuations near the critical point correctly. In particular, this 
method includes dynamical chiral symmetry breaking as the dominating fea- 
ture of QCD at low energies. The low-energy model lacks the gauge degrees 
of freedom, therefore it cannot be linked with lattice calculations for T > 160 
MeV. As a first step, we propose to extend the flow equations of the linear 
cx-model by a flow equation for free massless quarks and gluons at high scale. 
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Combining these two systems within a single renormalization group flow which 
integrates successively over both high and low energy excitations one obtains 
an adequate equation of state of strong interaction over a large range of tem- 
peratures and densities, as a comparison of our results with recent lattice 
simulations shows. 

The SU(2)l® SU(2)ji invariant linear cr-model with constituent quarks q and 
chiral mesons $ = (a, 7?) is an effective model for dynamical spontaneous chiral 
symmetry breaking at intermediate scales of < k < Am- It is given by the 
Euclidean Lagrangian density 

UoM = 0q + gq{v + irir^q + \{d^f + U{&) . 

The grand canonical partition function Z for the quantum theory has the path 
integral representation 

Z = Tre -/3(H-„N) (1) 

= jDqjDqjD® exp ^- J^dt E J d 3 x E {C LaM ~ /V?)) , 

where periodic and anti-periodic boundary conditions apply for bosons and 
fermions, respectively. A Gaussian approximation to the path integral followed 
by a Legendre transformation yields the 1-loop effective action 

H$ q q] =5-[<|> q q] -Trlog ( 6 -^M] + I Tr log ( l^Ml±) (2) 

Here, the boundary conditions of the functional integral appear in the mo- 
mentum traces and we neglect contributions from mixed quark meson loops. 
We consider the effective action T in a local potential approximation (LPA), 
which represents the lowest order in the derivative expansion and incorporates 
fermionic as well as bosonic contributions to the grand canonical potential 
density Q which is related to the pressure of the system by 

n(T,») = -p(T,n). (3) 



To derive RG flow equations, we use Schwinger's proper-time method to reg- 
ularize the respective logarithms. A straightforward partial evaluation of the 
trace in momentum space and a transformation of the expression to proper 



2 



time form 1 yields 



Q =Quv + I | ^ T g | ^ 47Vc7V/e - T ((, n+ ^) 2+p - + M|) (4) 



?1= — oo 



where the Matsubara frequencies take the values u n = 2nirT for bosons and 
u n = (2n+ l)nT for fermions, respectively. The effective masses are defined as 

M 2 =g 2 ® 2 , M 2 = 2^-, M 2 = 2^- + A^ 2 -^-. 

A general set of smooth cutoff functions f(rk 2 ), parametrized by a real positive 
number a fulfills the required regularization conditions, cf. ref. [5,6,9] 

In previous works [5,6,9], cutoff functions with integer a have been applied to 
various problems. However, in the case of a thermal system half-integer cutoff 
functions are better suited, since they effectively regularize only the spatial 
momentum integral. The relevant derivative of the cutoff function with respect 
to the IR cutoff scale 

k dk " r>TT) ( j 

involves a half-integer power of the proper time variable r, which combines 
with the half integer power from the integration over the spatial momentum 
coordinates. With a = 3/2 the proper time integral gives an analytic RG 
improved flow equation for the grand canonical potential Q 



u dn LaM (^,k) k 5 f 3/1 \ 1/1 ; 

(1 - n F (E q ) - n F (E q ))) -(6) 



2N c N f 



The appearing effective energies are defined by 

Ei = ^Jk 2 + M 2 , % G {tt, a, q} , 



1 The transformation of the apparently complex expression to proper time form is 
possible, since all imaginary parts cancel each other in the Matsubara sum. 
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and the occupation numbers have the usual form 

nB ^ = ^mzri> np W = e(E J /T + 1 . *f{e) = e(E J /T - i . 

Note the exceptionally simple form of this evolution equation. One can imme- 
diately read off the contributions from the vacuum and thermal dynamics of 
the respective particles appearing with their proper degeneracy factors. Espe- 
cially, compared to previous vacuum equations [9] energy denominators appear 
instead of the usual propagator terms, which is natural for a thermal equation 
because of the broken Poincare invariance. The analytic form of these equa- 
tions allows a numerical solution, especially in the case of finite temperature 
and density where the computational effort rises considerably. 
The thermodynamics and phase structure of the chiral system has been dis- 
cussed in [5-8]. One obtains a good description of the second order chiral 
phase transition at T c pa 160 MeV and also a qualitative picture of the first 
order phase transition at high densities. However, the model lacks the correct 
high-temperature behaviour since it does not contain the appropriate degrees 
of freedom in this regime. 

The fundamental degrees of freedom at high energies are quarks and gluons 
of QCD, described by the Euclidean Lagrangian density 

Cqcd = q(P + m c )q + \F^F^ , (7) 

where m c is the respective quark mass matrix. Due to asymptotic freedom 
quarks and gluons decouple at asymptotically high energies and the partition 
function takes a similar form as Eq. (1) with chiral mesons replaced by gluon 
fields. Lattice gauge simulations, cf. ref. [1,2], suggest that even at moderate 
temperatures T > 200 MeV the equation of state of the interacting system 
behaves almost like a free gas. Therefore in a first attempt to link the quark 
meson model with QCD we neglect the full gauge interactions and approxi- 
mate the high-energy part of QCD at scales k > Am by a free gas of massless 
quarks and gluons. The integration of the thermal fluctuations within a com- 
mon RG flow using the same smooth infrared cutoff function allows a smooth 
transition of the respective degrees of freedom. The matching depends on the 
chosen cutoff function. A smooth cutoff function has the same effect as a self- 
interaction of the gluons and provides an effective screening for the gluonic 
degrees of freedom. In order to match the two systems independently of the 
matching scale one would have to add the relevant quark-gluon dynamics. In 
this paper we show that, however, already the simple addition of the free quark 
gluon flow with the meson flow gives reasonable results for certain quantities 
as the equation of state. In further work we plan to apply the evolution of the 
gluon system [13-15,17,18] to reduce the dependence on the matching scale. 

We neglect the current quark masses for the two light quark flavors, then the 
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equation for the high-energy part has a very compact form in our simplified 
approximation: 



k 



QCD 



(k) k 4 



dk 



6n 2 
= k 4 



N g { l -+n B {k))-2N c N f {\ 



riF{k) — riF{k)) 



*e{s,<j,<j} 



Ni_ 
6ir 2 



+ Ui(k) 



(8) 



The sum runs over gluons, quarks and antiquarks including the appropriate 
occupation numbers with degeneracy factors and sign factors N g = 16, s g = +1 
for the bosons and N q = N q = 2N c N f =12, s q — Sq — — 1, for the fermions and 
antifermions, respectively. The potential of the total system is obtained by 
integrating the quark-gluon evolution equation from k = A^ to k = Am and 
the quark-meson evolution equation from k = Am to k = 0, i.e. the IR-cutoff 
parameter of the quark-gluon potential is identical to the UV-cutoff parameter 
of the quark-meson model. We can integrate the simplified QCD flow equation 
analytically and get 



&>qcd{Aoo, A m , T) —^2 s '> 



6tt 2 



k 4 



k 3 TLii (sinf(fc)) 



(9) 



-3k 2 T 2 Li 2 (sin°°(k)) - QkT 3 Li 3 (si<°(*0) - 6T 4 Li 4 (s i nf > (ife)) 



The polylogarithmic functions Li n (x) are defined by 

LiAx) — — m(l — x) , Li n (x) — [ — " 1 dz 

Jo z 

and the functions n°°(k) denote the classical Maxwell-Boltzmann distribution 
functions of the respective particles 

n™(k) = e~* , n™(k) = e"T^ , nf(k) = e" 4 ** . 

We can evaluate eq. (9) in the limit Aoo — > oo by making use of the behaviour 
of the polylogarithms in the vicinity of zero 

AT. ( A 4 A4 

fi Q CD(Aoo, Am, T) = ^ 8^ — f + -f - A M TLi 1 (s i n°° (Am)) (10) 
-3A 2 M T 2 Lt 2 (s t n°°(A M )) - QA M T 3 Lt 3 (s t n°°(A M )) - 6T 4 L2 4 (s J n° (A M ))) . 

Here we kept the divergent part proportional to A^, so the potential still 
needs to be regularized. In the limit of zero temperature the vacuum part can 
be identified, since it is only a function of Am and A^ 

AT. /A4 A4 \ 

^(Aoo, Am, T = 0) = £ ^ \-f - -fj . (11) 
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100 200 300 400 500 600 

T[MeV] 

Fig. 1. Different contributions to arising from the mesons and quarks (constituent 
quarks and free quarks ) (dashed line) and the gluons (dash-dotted line). The dotted 
line shows the scaled pressure function of the linear cr-model with two quark flavours 
which contains the cutoff Am- 



To obtain the potential QlctM of the chiral low-energy part, we integrate the 
corresponding flow equation from k = Am to k=0. As initial values we choose 
Am = 1.0 GeV, £/($ 2 ) = \m M <$> 2 + \\ M <$> A with \ M = 37, m M = 320 MeV 
and the quark-meson coupling g = 3.45. These initial parameters are adjusted 
to get the value of the pion decay constant < a >= f n = 88 MeV at k=0 in 
the chiral limit. In principle, the starting parameters should result from the 
integration over the QCD degrees of freedom. This may be possible including 
quark interactions from gluons at higher scales. 

The full potential reads: 

*WAoo, k = 0, T) = fi QC D(Aoo, A M , T) + n LaM (A M , k — 0,T). (12) 

The finite-temperature part relevant for the pressure does not depend on the 
vacuum energy. Therefore in the regularized potential the UV-divergent part 
drops out as it must 

p(T) = -n total (A oo ,k = 0,T) + n total (A oo ,k = 0,T = 0). (13) 

In the limit of infinite temperature T — > oo the ratio of p/T 4 yields the Stefan- 
Boltzmann limit for a massless quark-gluon gas, because the region k < A M 
has vanishing weight 



y4 



N 9 + l(N q + N 9 j) ^ . (14) 
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Fig. 2. .The pion decay constant as a function of temperature for vanishing chemical 
potential. 



For small temperatures T — > we find a massless gas of pions with 

P(T) 



rp4 



3^. 
90 



(15) 



In fig. 1 we show the different contributions to p/T 4 arising from the integra- 
tion over the flow of the mesons and all quark contributions (including both 
the constituent quarks with cutoff Am and free quarks without cutoff) (dashed 
line) and the integration over the gluons (dash-dotted line). Furthermore we 
plot the scaled pressure of the linear a-model with cutoff A M (dotted line) 
described by eq. (6). 

The comparison indicates clearly that the quark-meson model with cutoff 
A M alone becomes unreliable for a description of the equation of state at 
temperatures T > 150 MeV. It definitely has to be linked with quark and 
gluon degrees of freedom at higher temperatures. The linear sigma model is 
adequate at low temperature where massless pions dominate the pressure and 
the quarks are massive due to spontaneous breaking of the chiral symmetry. 
Gluons contribute to the pressure only at rather high temperatures since the 
smooth infrared cutoff around the scale A M introduces a gluon mass and gives 
deviations between our model and the free gluon gas. The flow equation with 
the successive integration over different degrees of freedom avoids a first-order 
transition which typically arises in simplified models with two components like 
the bag model combined with a pion gas [16]. 

The critical behaviour of the system is described by the pion decay constant 
which is shown as a function of temperature for vanishing chemical potential 
in fig. 2. The figure shows that the pion decay constant goes to zero for T c = 
154 MeV continuously. For temperatures below T c the quarks are massive 
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Fig. 3. The scaled pressure density multiplied by a factor of three compared to 
the scaled energy density as a function of temperature. The upper and lower Ste- 
fan-Boltzmann limits corresponding to a free quark-gluon gas and a pion gas are 
drawn as dashed lines. 

and the pions are massless, so the system is in the chiral broken phase. For 
temperatures above T c the system is in the chiral symmetric phase, where 
the quarks are massless and the pions are massive. So one obtains a second 
order chiral phase transition at T = T c , which is about 20 MeV smaller than 
the result of lattice calculations [2]. Note that the pion decay constant as a 
function of temperature for the chiral low energy system does not differ from 
the function of the total system, because the pion decay constant is zero in 
the high-energy part of our model. 

In fig. 3 we compare the curve of the scaled energy ^ with three times the 
pressure ^ as a function of temperature. Both quantities || gradually 
increase across the second order phase transition. The scaled pressure density 
reaches about 95% of the Stefan-Boltzmann limit at T pa 4 T c . Respective 
lattice calculations [2] reach about 85% of the Stefan-Boltzmann limit at T pa 
4 T c . The energy density becomes more rapidly asymptotic, cf. fig. 3. Thereby, 
the "interaction measure" e— 3p has a peak at T pa 1.27 T c pa 196 MeV, whereas 
the respective maximum of lattice calculations of (2+l)-flavour QCD [3] is 
located at approximately T pa 1.1 T c pa 189 MeV. The "interaction measure" 
e — 3p has more subtle features at finite chemical potential. For instance we 
find a low temperature maximum in e — 3p at finite chemical potential which 
does not agree with finite density lattice calculations. This low temperature 
maximum arises from the fact that we have constituent quarks in our model 
and no baryons for temperatures below T c . Lattice calculations of the same 
quantity go rapidly to zero for temperatures below T c which may be an effect 
of confinement and/or of the limited grid sizes and the high masses of the 
nucleons. 



8 



100 200 300 400 500 600 

T[MeV] 

Fig. 4. The scaled pressure difference = p ^ T '^~p^ T '^~ ^ [ s shown as a funtion of 
temperature for /j>b = 100,210,330,410,530 MeV (from bottom to top). 

To study the effects of a finite baryon chemical potential on the pressure it is 
useful to consider the difference Ap of the pressure of the baryonic system to 
the nonbaryonic system at the same temperature: 

Ap=(p(T,fi)-p(T,fi = 0)). (16) 

Small chemical potentials compared to the baryon mass suppress the baryonic 
pressure at low temperatures. Our calculation shows a suppression of Ap/T 4 
in agreement with the low-temperature expansion for a gas of nonrelativistic 
and non-interacting constituent quarks with masses m > /i = /is/3, which 
reads 

Ap g / m s 



T 4 2V2nl VT 



-e"f . (17) 



Here g denotes the degeneracy factor of the fermions. In fig. 4 one can see the 
scaled pressure j£ of the quarks for different chemical potentials //# = 100 
MeV to /iB — 530 MeV. We have chosen these values identical to the ones in 
the lattice calculation in ref. [3], where the equation of state includes heavy 
strange quarks. At high temperatures the scaled pressure difference decreases 
according to 

^ = (-!*- + -J* V (18) 

T 4 } \24T 2 48tt 2 T 4 J y ' 

This behaviour is also reproduced by our calculation, cf. fig. 4. 

Linking the quark-meson model with QCD degrees of freedom we find that the 
low-energy sigma model can be extended to high temperature T > 150 MeV. 
The combination of different degrees of freedom within the renormalization 
group flow equations opens the window for a consistent description of QCD 
both at low and high scales which is needed in the equation of state. Already 
the simplified model we constructed shows good agreement with lattice simula- 
tions. As indicated above, the low-momentum behaviour of the model has to be 
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improved to get better agreement with the physics for temperatures below T c . 
Furthermore the matching between the quark-gluon and quark-meson model 
is too simple and should be extended by including the relevant quark-gluon 
dynamics on the high-energy side. Both the effective four-fermion coupling 
and gluon condensation depend on the matching scale. A scale independent 
transition may be achieved if one can connect the chiral symmetry dynamics 
with gluon condensation. 
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